Semiclassical magneto-transport in graphene n-p junctions 
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We provide a semiclassical description of the electronic transport through graphene n-p junctions 
in the quantum Hall regime. This framework is known to experimentally exhibit conductance 
plateaus whose origin is still not fully understood. In the magnetic regime {E < vpB), we show 
the conductance of excited states is essentially zero, while that of the ground state depends on the 
boundary conditions considered at the edge of the sample. In the electric regime (E > vfB), for a 
step-like electrostatic potential (abrupt on the scale of the magnetic length), we derive a semiclassical 
approximation for the conductance in terms of the various snake-like trajectories at the interface 
of the junction. For a symmetric configuration, the general result can be recovered using a simple 
scattering approach, providing a transparent analysis of the problem under study. We thoroughly 
discuss the semiclassical predicted behavior for the conductance and conclude that any approach 
using fully phase-coherent electrons will hardly account for the experimentally observed plateaus. 

PACS numbers: 73.22.Pr, 73.43. Jn, 03.65.Sq, 73. 23. Ad 



I. INTRODUCTION 

Graphene, a two-dimensional material made of a 
mono-layer of carbon atoms arranged in a hexagonal lat- 
tice, has been studied theoretically for some time [IHl], 
but means to isolate and manipulate it experimentally 
have only been developed a few years ago |S]. Since 
then, the exotic properties of this material have aroused 
considerable interest in the condensed matter physics, 
chemistry, and material science communities [6l|7]. These 
properties mostly originate from graphene's peculiar low- 
energy band structure, both gapless and linear. Con- 
duction and valence bands touch each other at the cor- 
ners of the (hexagonal) first Brillouin zone, two of which 
{K,K') are inequivalent. At low energies, the dispersion 
relation vanishes at these Dirac points and the effective 
Hamiltonian for electrons in graphene is that of a mass- 
less pseudo-relativistic particle travelling at an effective 
"speed of light" vp- Here, the internal degree of freedom 
associated with graphene's two inequivalent sublattices 
plays the role of the spin (the true spin providing merely 
a degeneracy factor as long as spin-orbit interaction or 
magnetic impurities are neglected, as we shall do here). 

Two of the most prominent features that distinguish 
graphene from conventional two-dimensional electron 
gases (2DEG) are the occurrence of Klein tunneling [51- 
[TT| . i.e. the absence of backscattering of charge carriers 
when approaching potential barriers under normal inci- 
dence, and that of an anomalous quantum Hall effect 
(QHE) |12|, I13j stemming from the peculiar quantization 
of the Landau levels in this material. 

A nice experiment probing the physics that emerges 
from the interplay of Klein tunneling and the anomalous 
QHE was carried out in 2007 by Williams and collab- 
orators [H]. Using a global back gate and a local top 
gate, they were able to apply gate voltages of different 



signs in two regions of a graphene sheet, creating a n-p 
junction at the regions interface, an ideal setup to study 
Klein tunneling. Using a strong perpendicular magnetic 
field they performed quantum Hall transport measure- 
ments across this graphene n-p junction. Quite remark- 
ably, for some combinations of n and p filling factors, the 
experimental data exhibit conductance plateaus at the 
fractions 1/2 and 3/4 of the quantum of conductance of 
spin-degenerate systems, Gq = 2e^ /h, at odds with the 
usual graphene quantum Hall series. 

These observations were confirmed by other experi- 
ments p!5Hl7] and extended to the setup of single-layer 
n-p- n junctions (TS] and bilayer n-p-n junctions [19j. The 
observed plateaus in n-p junctions can be cast in terms 
of the QHE edge channel picture, within good accuracy, 
by the expression 

where iV„ and Np stand for the number of incom- 
ing/outgoing edge modes, which depend on the gate volt- 
ages applied to the n and p regions, and consequently on 
their filling factors. 

Very early, Abanin and Levitov [20] proposed an in- 
terpretation for Eq. ([I]) in terms of a "quantum chaos 
hypothesis" . The key element is that since n and p edge 
channels in the quantum Hall regime have opposite chi- 
rality (direction of propagation) , they interfere with each 
other at the junction interface. The "quantum chaos hy- 
pothesis" amounts to assume that this interference ef- 
fect gives rise to a mode-mixing mechanism that is suf- 
ficiently strong that an incoming channel at the inter- 
face will have equally likely probabilities of leaving the 
interface in any of the outgoing channels, leading to a 
current partition consistent with Eq. ([I]). This picture is 
further motivated by the observation that the values of 
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the conductance plateaus Gnp coincide with the average 
conductance (G) , predicted by the random matrix theory 
(RMT) for chaotic systems [^[^ . 

Albeit very appealing, a complete interpretation of the 
experimental results by means of the quantum chaos hy- 
pothesis is still problematic. As already pointed out in 
Ref. 20], RMT gives an average (G) = G„p, but also pre- 
dicts universal conductance fluctuations (UCF) of the or- 
der of e^//i. Those are definitely not currently observed 
by experiments. UCF are ubiquitous in both quantum 
chaotic and diffusive systems [H Hi]. They 

have been studied experimentally as well as theoretically 
for ordinary graphene flakes (see, for instance, Ref. |25j 
for a review). In the quantum Hall regime, numeri- 
cal investigations of graphene n-p junctions [27] and 
of graphene-superconductor junctions [28] showed that 
modifications of the disorder strength may cause small 
changes in the magnitude of the UCF but do not im- 
ply qualitative modifications of the transmission. These 
observations thus support the conclusion that in the ab- 
sence of decoherence, the role of which still has to be es- 
tablished in the problem we are investigating, disorder by 
itself cannot explain the conductance plateaus observed 
experimentally in quantum Hall n-p junctions. 

Further insight is gained from another recent experi- 
mental work investigating the fingerprints of Klein tun- 
neling in a graphene p-n-p junction, this time in the 
presence of a weak magnetic field [TT]. In this case, the 
authors observe beautiful quantum interference patterns 
quantitatively explained by a semiclassical approxima- 
tion in terms of a set of simple transmitted and reflected 
paths across the junction. It is noteworthy that, although 
the sample quality of all cited experiments is similar, in 
Ref. [TT] there is no indication that disorder plays a sig- 
nificant role. Obviously, one should be cautious before 
drawing conclusions from this observation, since by go- 
ing from small to large magnetic fields one dramatically 
changes the physics of the problem. Nonetheless, it cer- 
tainly justifies the study of the transmission properties 
of clean n-p junctions in the quantum Hall regime. 

An even stronger motivation to consider a model of 
clean ballistic graphene is that, at the classical level, 
the snake-like [50] trajectories of electrons on a graphene 
step-like symmetric n-p junction have equal refiection 
and transmission probabilities (see beginning of section 
IV), leading to current partition consistent with the ex- 
periments. 

The goal of this paper is to use a semiclassical analy- 
sis to investigate whether the classical current partition 
behavior is transposed into the "quantum chaos hypothe- 
sis" in the quantum coherent transport limit. Extending 
over previous results already presented in [29], we shall 
see that the answer is actually negative, and that clean 
graphene n-p junctions in a high magnetic field provide 
an example where the classical behavior is not a good 
approximation to the semiclassical (or quantum) behav- 
ior. We find that, in the full quantum coherent transport 
regime, the transmission fluctuations due to quantum in- 



terference between different electronic paths across the 
junction are large and non-universal. The relation be- 
tween these results and the experimental ones is discussed 
in sections \VJ\ and IVlII 

Our analysis is based on the Landauer conductance 
formula 



G — Go T„ 



(2) 



where T„_q, is the transmission coefficient, n stands for the 
channel index, and a for the valley index that specifies 
in which of the valleys K {a = or K' {a = —1) the 
edge channel is polarized [HI] . We evaluate the transmis- 
sion coefficients via the Green's functions formalism in- 
troduced by Fisher and Lee [30] and later generalized by 
Baranger and Stone |31J (mainly to include the presence 
of a magnetic field). Here we use semiclassical Green's 
functions obtained for graphene by two of us in a previous 
publication [32]. The semiclassical approximation allows 
us to express G analytically and can be used to model 
systems of realistic sizes, which are hardly reachable by 
numerical techniques that use an atomistic basis. 

The description of a graphene n-p junction depends 
sensitively on the relative strength of the magnetic and 
electric fields at the interface between the electron-doped 
and "hole" -doped regions. This is best understood by 
recalling the insight gained in graphene electronic trans- 
port in presence of crossed electric (in plane) and mag- 
netic (perpendicular) fields [SHI 131] : the electromagnetic 
Lorentz invariant quantities Ci = E • B and £2 — 
— w|nB^ define two distinct transport regimes char- 
acterized by the parameter 



/3 
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If /3 < 1 the electric field will always dominate the mag- 
netic field, as it is possible to apply a Lorentz boost and 
go to a reference frame where the magnetic field is made 
zero. This is referred to as the electric regime. In con- 
trast, if /3 > 1 the electric field can be eliminated by go- 
ing to a proper reference frame. This regime is referred 
to as the magnetic one. The Lorentz invariance stated 
above relies on translational invariance, and hence is not 
expected to hold near the edges of a graphene ribbon. 
Nonetheless, the classification in terms of electric and 
magnetic regimes still makes sense for confined systems, 
as will be shown in the following. 

The remainder of the paper is organized as follows. 
We begin presenting our model and the key elements of 
the semiclassical analysis. In section |III[ we then dis- 
cuss the physics of the magnetic regime (in particular 
the limit /3 ^ 1). Within this regime (for which no quan- 
tum chaos hypothesis has been evoked), we show that a 
classical description based on the concept of adiabatic 
invariance suffices to explain most of the junction's fea- 
tures. We obtain that the conductance is essentially zero 
for all modes except for the lowest one, whose transmis- 
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sion depends on the boundary conditions considered at 
the edge of the ribbon [35l |36] . 

The electric regime, which we beheve to be that of the 
experiments [rH - fTBl [T5] , is addressed in sections IV and 
\V\ Those constitute the main technical development of 
this paper. At their end, we present a summary and a 
discussion of the main results they contain. In section 
|IV| the electric potential is taken to be symmetric, such 
that the (absolute) doping is the same at the n and p 
regions. This allows for a simple and transparent anal- 
ysis using a semiclassical scattering matrix approach. A 
general treatment, beyond the symmetric case, of the 
electric regime is presented in section |Vj It requires de- 
termining the exhaustive list of trajectories involved in 
the dynamics at the interface, and making use of the 
semiclassical approximation to the single particle Green's 
function in graphene |32j to compute the Fisher-Lee for- 
mulae. We present generalized expressions of the latter, 
which take into account the internal pseudo-spin struc- 
ture of charge carriers in graphene. Some qualitative fea- 
tures of the transmission in the general case are similar 
to the ones found in the simpler symmetric configuration. 
In particular, in both cases the conductance differs sig- 
nificantly from the classical one. A qualitative analysis of 
the robustness of our semiclassical findings with respect 
to modifications of the model such as boundary condi- 
tions, disorder and steepness of the potential barrier is 
conducted in section 



VI Our conclusions are presented 



in section rVTT 



II. THE MODEL 

We consider a graphene sheet of finite width W in the 
y direction, and connected in the x direction to infinite 
leads on both sides. We assume the n-p junction is ob- 
tained by applying an electrostatic potential in the plane 
of the graphene sheet, separating it in three distinct re- 
gions: 

{Vi , X <0 (region 1) 
Vi + eEx , 0<x<L (region J) , (4) 
V2 = Vi + A , x > L (region 2) 

with Vi < 0, V2 > (we assume the chemical potential 
fixed at 0), and A — eEL the height of the potential 
step. L governs the steepness of the potential step and 
E is the intensity of the associated electric field E = Eic, 
where x is the unit vector in the x direction. Applying 
a strong perpendicular magnetic field B = Bz, suc h that 
the quantum Hall regime Ib<^W {Ib = ^h/{eB) is the 
magnetic length) is reached, leads to the classical picture 
of skipping-orbit motion at the edges of regions 1 and 2, 
with a direction of propagation depending on the edge 
and the type of charge carriers. The electron dynamics 
in region S and in its neighborhood is what interests us 
here. 



To describe the quantum dynamics of the electrons 
in the graphene sheet, we use the representation \I/ — 
{'4'At4'Bt~'4''bi'4''a)'^ ^^'^ wave function. Here, as 
usual, {A, B) denote the real-space graphene sublattices, 
{iPaj^b)'^ the wave-function in valley K and ("0^, '0^)'^ 
that in valley K' . In such a representation the low-energy 
Hamiltonian of graphene electrons is isotropic in valley 
space and takes the form 

y. = vpTo <Si{Il-d) + U (r)To (To + m{r)vpTz ® CTz , (5) 

where C/(r) is the electric potential generating E, the 
magnetic field B derives from the vector potential A(r) 
via the Peierls substitution 11 = p -I- eA(r), and a mass 
term m(r) has been included for completeness. In Eq. ([s]) 
T and (7 are Pauli matrices respectively associated to val- 
ley space and sublattice space (with the usual convention 
that To = (To = 1), and the scalar product is restricted 
to the x-y graphene plane. 

Three types of boundary conditions will be discussed: 
most of the time we shall assume infinite mass confine- 
ment, but also in some circumstances zigzag and arm- 
chair. All of these can be incorporated in a matrix 
representation introduced by Akhmerov and Beenakker 

[371 133, 



M = ■ f) {fi ■ a) , 



(6) 



where n is the normal to the outward pointing unit vec- 
tor at the boundary and ? is the polarization of the edge 
state in valley space. Each boundary condition is asso- 
ciated to a given matrix A4, which actually amounts to 
a given polarization v. For instance, the case of infinite 
mass confinement corresponds to imposing the condition 
'^b/'^a — ^'b/'^'a — *exp(i7) on the boundary with 
7 — cos~^(n-y), as shown by Berry and Mondragon |39| . 
This is obtained by choosing i7 = z in Eq. ^. Similar 
parameterizations are obtained for zigzag and armchair 
lattice terminations (see Refs. [571 138) ). 

At the leads, that is, within regions 1 and 2 of Eq. Q, a 
propagating mode (n, a) in the quantum Hall regime can 
be built on the set of skipping trajectories emerging from 
the boundary with an angle 0„ ^ such that the transverse 
action is quantized, according to the generalized EBK 
formula [40j 



dr 



-6 



27rn 



(7) 



In this quantization condition, the closed path on which 
the integration is performed can be taken as the one 
shown on Fig. [I] (bearing in mind that on r2 the path 
does not follow a trajectory, and thus the momentum is 
not parallel to the shown direction), /i is a Maslov index 
counting the caustics traversed by the closed path and 
(/)fcc is the phase acquired by the electron upon bouncing 
on the boundary. For graphene, one needs in addition to 
take into account a semiclassical phase ^sc which, in the 



4 




FIG. 1: (Color online) Illustration of a closed path F — riUr2 
along which the action integral in Eq. ([7| is computed. The 
phase of the mode wave-function must remain unchanged (up 
to an integer multiple of 2n), leading to the quantization of 
the edge angle 0„^a- 



absence of a mass term in the Hamiltonian ^ coincides 
with a Berry phase p2l, and is just given by half the 
rotation angle of the momentum vector (see the discus- 
sion at the beginning of section |v| . Because inside the 
graphene ribbon we only consider a massless Hamilto- 
nian, in what follows we refer to the semiclassical phase 
as a Berry phase, for simplicity. 

For infinite mass boundary conditions, the phase ac- 
quired upon reflection on the boundary (phc = (f'moa 
Eq. ([7]) can be easily shown to be (t)m^ = tt (Dirich- 
let) in valley K and 4>m^ = (Neumann) in valley K' . 
The path of integration F traverses one caustic (so that 
/i = 1) and = — . Wrapping everything up, the 
quantization condition can be written as 



sin 29,, 



= 27rS(r 



4' 



with 



Ri — 



eBvp 



(8) 



(9) 



the cyclotron radius at the electron {i = 1) or hole {i = 2) 
side. Note the quantization condition ([s]) with n ~ can 
be fulfilled for only one of the two valleys {a = —1, for 
which Neumann boundary conditions apply). Introduc- 
ing the function 



/ : e -^0-sm{29)/2 
and noting that eBR^/h = 2vi, where 



(10) 



(11) 



is the filling factor at the electron [i = 1) or hole {i ~ 2) 
region, we recast Eq. (Isl) as 



n 4 



(12) 



This form of the quantization condition makes explicit 
that, besides the indices {n,a), the quantized angle 9n^a 
only depends on the filling factor Vi. Quantization con- 
ditions for zigzag and armchair boundary conditions can 
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FIG. 2: (Color online) Skipping orbit motion for classical 
electrons coming from the upper left edge with an angle 6 = 
7r/4, followed by bubbling. Different panels correspond to 
various values of j3 in the magnetic regime. The black dot- 
dashed lines indicate where the junction region begins and 
the red dotted lines indicate where bubbling is expected to 
occur according to the semiclassical theory. The adiabatic 
transition from the edge to the bulk is apparent in the upper 
panels, while it clearly breaks down in the lower ones even 
though classical bubbling still occurs. 



be derived in the same way |3Tj . Note finally that as the 
angle 9n.a — ^ ti", which corresponds to a transition to a 
Landau level in the bulk (for which a can be interpreted 
as switching to 0) a uniform approximation should be 
used. This uniform approximation is described for the 
scalar (Schrodinger) case in Ref. [42) . 



III. MAGNETIC REGIME 

In this section, we study the transport in the magnetic 
regime in the limit of /3 3> 1. This corresponds to the sit- 
uation where the electric potential varies adiabatically on 
the scale of the cyclotron radius. One can then use that 
the action ([s]) of the electron along the edge is an adia- 
batic invariant, and is thus conserved while propagating 
in the x direction. Defining an effective local cyclotron 
radius R{x) = \V{x)\/{eBvF), conservation of the ac- 
tion implies that the decrease in R{x) must be compen- 
sated by an increase in the effective angle 9{x) along the 
edge, since the function f{9) in Eq. (10) is monotonous 
on the interval [0,7r]. When this angle reaches the value 
TT, a transition takes place from the edge skipping-orbit 



motion to the transverse bulk cyclotronic motion with a 
vertical drift caused by the electric field. We will hence- 
forth refer to this transition by the visually appropriate 
term of "bubbling" . Once an electron has bubbled from 
the edge, it traverses the graphene ribbon transversally 
until it reaches the opposite edge, and then propagates 
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Vj = 1 V, = 




Vj/Eg W/a 



FIG. 3: (Color online) Exact (quantum mechanical) conduc- 
tance G across the junction as a function of the potential V2 
(in units of es = hvp/ls) in the right side of the junction, 
with a fixed value —V\/tB = 1-82 of the potential in the left 
side of the junction (which in terms of the filling factor is 
1^1 — 1.65). Solid line (black online): "infinite" mass bound- 
ary condition (note the mass actually used in the numeri- 
cal simulation is m = 0.5t, with t the hopping parameter); 
Dashed (red online): armchair with W/a = 80 (a is the dis- 
tance between nearest carbon atoms of the same sublattice); 
Dotted with rhomb symbols (green online): armchair with 
W/a = 81; Dotted with triangle symbols (blue online): arm- 
chair with W/a — 82. The vertical dot-dashed line indicates 
values of V2 for which the filling ratio V2 in region 2 is an 
integer. When V2 < the system is effectively a n-n junc- 
tion, and the conductance is given by the number of channels 
in region 2, namely 2[i'2] + 1 since min([i/i], [1^2]) ~ [1^2] here 
(see discussion in main text). For positive V2, the conduc- 
tance across the n-p junction is governed by the transmission 
of the lowest channel and, as expected according to Eq. ( 14 1, 



is given by G = Go for armchair boundary conditions with 
W/a = 2 (modulo 3), G = Go/4 for other armchair widths 
and G — for the infinite mass confinement case. The plot 
was obtained for /3 = 5, with a system of length L/a = 60 
and a magnetic length Is/a = 6.8. 



back into the lead that it came from. This scenario is 
illustrated in Fig. [2] 

In the case of a n-n junction, we see that the trans- 
mission probability of an electron can only be or 1. 
Quantitatively, the bubbling abscissa of an electron is de- 
termined by the relation Rff{0n,a) — R'^{xbub)T^, which 
leads to the value Xbub — (3Ri{l — \/ f{Sn,a)/T^)- A given 
mode (n, a) bubbles if Xbut < L, which can be equiva- 
lently expressed through the bubbling condition 



V2_ 

\Vi\ 



< 



f{On 



(13) 



FIG. 4: (Color online) Exact (quantum mechanical) conduc- 
tance G across the junction as a function of the width W (in 
units of a) for different values of the parameter j3 governing 
the transition between the magnetic (/3 ^ 1) and the electric 
(/3 <^ 1) regimes. The system considered was confined by an 
"infinite" mass m/t = 0.5 and the filling factors on both sides 
of the junction were taken to be — 1.65 and 1^2 = 0.10 
(hole doped). The length of the system and the value of the 
magnetic field are the same as in the previous figure. We ob- 
serve that the adiabatic result G ~ expected from Eq. ( |14[ | is 
observed for /3 = 5, but that /3 = 2 already shows some oscilla- 
tions in the conductance. For /3 < 1, the featured oscillations 
resemble those predicted by our semiclassical calculations in 
sections IIVI and IVl 



G = Gomm{Ni,N2) [20], with N, = 2[iy,] + 1 (with [•] 
the integer part), as expected in the context of quantum 
adiabatic transport [33]. Note that the quantum adia- 
baticity criterion, namely that the electrostatic potential 
must vary less than the inter-Landau level spacing on the 
scale of the magnetic length, is actually more robust that 
the classical one. 

For a n-p junction, transmission is zero for all modes, 
except possibly for the zero-energy mode for which the 
semiclassical reasoning above cannot be applied. The 
latter must be treated separately. Hence, a n-p junction 
in the adiabatic limit (3^1 has a conductance G = 
GqTq, with To the transmission probability of the n — 
mode. 

In the absence of inter-valley scattering, the transmis- 
sion coefficient Tq of a single channel through a n-p junc- 
tion in the quantum Hall regime depends only on the 
valley-polarization of the edge states [35] [3S] , namely 



Tn = 



1 — cos ( 



(14) 



(which of course is always realized for n-p junctions). 
This leads naturally for a n-n junction to the result that 



where the angle (f>i, is the one separating the valley- 
polarization vectors of top and bottom edge states on 
the Bloch sphere. For an armchair ribbon, this leads to 
plateaus in the conductance at Go or Go/4 depending 




FIG. 5: (Color online) Sketch of the geometry of a parallel 
junction. Due to the smooth transition from the edge of the 
ribbon to the n-p step-like junction, the quantization of the 
incoming channels is maintained. 



on the width of the ribbon [35]. For a zigzag ribbon, it 



was reahzed that formula ( 14 ) cannot be applied since 



a potential barrier, no matter how smooth it is, causes 
inter- valley scattering |36j . However, very similarly to 
the armchair case, it turns out that the transmission de- 
pends on the width of the ribbon and can be either zero 
or one &2\ . The case of infinite mass confinement can be 
treated similarly to the armchair one, yielding zero trans- 
mission just like for the higher modes. These results for 
the armchair as well as infinite mass confinement cases 
are illustrated on Fig. [3) The perfect refiection for in- 
finite mass confinement is a trivial consequence of the 
conservation of the valley-polarization of the edge state. 
The expected result (in the limit /3 ^ 1) for the conduc- 
tance of a graphene n-p junction laterally confined by an 
infinite mass is thus zero, up to exponentially small tun- 
neling contributions. As illustrated in Fig.[4j this limit is 
already reached from this point of view for /3 ~ 5, when 
a /? ~ 2 is already in the transition towards the elec- 
tric regime where some oscillations in transmission (and 
thus the conductance) are already visible. The data pre- 
sented in Figs.|3)|4]were obtained with a recursive Green's 
function technique, using the numerical software KNIT 
developed by Kazymyrenko and Waintal [44 . 



FIG. 6: (Color online) Typical trajectory along a parallel 
junction (rotated by an angle 7r/2 with respect to Fig. [5| 
between initial x = x' and final x — x" Poincare sections. 
Excursions cover a length Li on the electron side and L2 on 
the hole side. Angles with respect to the interface axis are 
shown. 



that quantization of the edge modes is maintained when 
arriving on the n-p interface. This amounts to taking 
the junction "parallel" to the edge of the ribbon instead 
of perpendicular to it. The main reason for making this 
choice is of course that the parallel junction is a simpler 
problem to tackle analytically. However, as discussed 
in more details in [45] . it can be shown that except for 
diffractive-like contributions at the edge-junction corner, 
going from a geometry for which the junction is perpen- 
dicular to the edges of the ribbon to one where it is paral- 
lel mainly amounts to applying a unitary transformation 
to the mode basis, under which the Landauer-Biittiker 
formula ([2| is invariant. As in any case a completely re- 
alistic description of the dynamics in the corners would 
depend on many details not included in our model (and 
probably unknown), the parallel junction is presumably 
as close (or as far) from a perfectly realistic description 
of the junction than a perpendicular one. 

We start by introducing a few notations for trajectories 
such as the one illustrated on Fig. [6] Let us denote by 
6*1 and 62 the angles between the x-axis and the vector 
n = p -f eA(r) (which is parallel to the velocity in the 
n region, but antiparallel to it in the p region) when the 
trajectory emerges from the edge or from the junction in 
respectively the n and p side. The two angles are related 
through the Snell-Descartes law 



|yi|cOs6li = 1^2 COS 6*2 



(15) 



IV. ELECTRIC REGIME: SYMMETRIC CASE 

We switch now to the electric regime. Two important 
simplifications are made in the next couple of sections. 
First, we consider the electric potential to be step-like 
on the scale of the magnetic length (but not on that of 
the carbon lattice, so as to avoid inter- valley scattering), 
placing ourselves in the opposite limit (/3 <C 1) as that 
of section [hi] This makes it possible to neglect the mag- 
netic field during the interaction with the barrier. The 
second simplifying hypothesis has to do with geometry. 
Compared to the experimental setup, we will consider a 
graphene ribbon where the transition from the edge to 
the step-like junction is very smooth (as in Fig. pi), such 



which expresses conservation of the momentum in the x 
direction. During the time interval between two consec- 
utive bounces on the edge boundary or at the junction 
interface (we call this an "excursion"), an electron cov- 
ers a distance Li = 2i?i sin^i on the n side (respectively 
L2 = 2i?2 sin02 on the p side) in the longitudinal x di- 
rection. 

Klein tunneling or refiection at the junction interface is 
given by the same probability amplitude as in the absence 
of magnetic field, namely [8] 

fi = e'^^ri for refiection 1 — > 1 (16) 

ti = e^(ei+e2)/2^^ f^^. transmission 1^2, (17) 
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with 



ti 



cos((gi + 02)/2) 
'cos((^i-02)/2) 
sin 6i 



(18) 
(19) 



For incident waves from the p side, the corresponding 
expressions for (^2,^2) and (^2,^2) a-re obtained by ex- 
changing the roles played by 9i and 09. The phases 
of the factors e^'^i and e'C'i+^^'Z^ in ([16) and ^\ can 
be interpreted as Berry phases as they correspond to 
half the pseudo-momentum (which as already mentioned 
is, for holes, antiparallel to that of the velocity) ro- 
tation during the scattering on the junction interface 
(01 = -(-01 - 0i)/2 and (01 + 02)/2 = -(-02 - 0i)/2). 
We have therefore distinguished them from the "genuine" 
reflection and transmission coefficients ri and ti given by 



( 18 1 and ( 19 ) 



The transmission probability through the interface is 
deduced from the quantum amplitude by taking into ac- 
count the flux normal to the barrier: 



„ sin 09 , ,9 sin 01 , 

sm 01 sm 02 cos 



sin 01 sin 02 
M(^?i-^2)/2) 



(20) 



The rest of this section will be devoted to the symmet- 
ric case V2 = —Vi. This leads via Eq. (|9]) to i?i = R2 and 
via Eq. (15) to 02 = 0i, and hence L2 = Li. Excursions 



on both sides of the interface cover the same distance, 
implying that the reflected and transmitted waves of a 
scattering charge carrier meet at equidistant "vertices" 
(see Fig. [7]). 

Let us start by quickly discussing what is expected 
classically for this configuration. A classical incident 
electron has a probability T = sin^ of being trans- 
mitted through the (symmetric) potential step. Calling 
Xi = (cijhi)'^ the vector composed of the probabilities 
Ci and hi for the incoming particle to emerge at vertex i 
on the n or p sides of the junction (e^ + hi — I) we have 



Xi+i — 



l-T T 
T l-T 



X,, 



(21) 



The matrix in Eq. ( 21 1 can be diagonalized and has eigen- 



values 1 and A = 1 — 2T, which leads to 



X 



N 



1 + A^ 1-A^ 
1-A^ 1 + A" 



(22) 



Since |A| < 1, the asymptotic behavior is as expected 
Xoo = (1/2,1/2)'^ and Eq. ^ tells us this limit 



is reached exponentially quickly with the number of 
bounces on the interface. An incoming classical electron 
has thus equal chances of being reflected or transmitted, 
provided the interface between regions 1 and 2 is long 
enough. Let us now show that this is no longer the case 
for a quantum particle. 

Semiclassically, one needs now to propagate the ampli- 
tudes z*^*^^ and z^'*-' on the electron and hole sides from 




FIG. 7; (Color online) In the symmetric case, the length L 
of excursions on both sides of the junction are equal, making 
reflected and transmitted waves meet at equidistant vertices 
along the interface. N is the largest number of excursions for 
a given interface length W . 



one vertex to the other. Noting Zi — {zf \ zf^^Y" these 
amplitudes at vertex z, this propagation can be obtained 
as Zi-^i = SZi, where the scattering matrix can be writ- 
ten as a product S — VTZ with 



V = 



f Si— if/Ji+i^i 





eA'S'2 







(23) 



describing the propagation in the n and p regions and 



7^ 



h 7=2 



(24) 



the transmission or reflection taking place at the in- 
terface. The matrix V implies mainly a multiplica- 
tion by a phase, which includes the action integral Si 
along the classical trajectory, the Maslov phase /i^ as- 
sociated with the traversal of caustics, and the Berry 
phase ^i associated with (half) the rotation of the pseudo- 
momentum vector n. One obtains for these various quan- 
tities 5*1 — eBR^f{dn,a), with the function / defined in 
Eq. (10), 5*2 = 5*1 — 2'!Thh' where v = vi = 1^2 is the 



filling factor (11) in both the electron and hole regions, 
^1 = ^2 = —0n,a and fj,i = —pL2 — 1 (the Maslov index 
accounts for a single caustic and is counted negatively 
on the hole side since velocity and momentum are op- 
posite there). For the symmetric junctions we consider 
here, we furthermore have f 1 = 7=2 = — e*^"'° cos0„ „ and 



so that finally 



-ie^'^^ cos 0„ a 



e"''' sin 



— e *'^'^sin0„,a ie 



cos On, a 



(25) 

For a given channel (n,Q;), the global phase factor in the 
scattering matrix expression is irrelevant and will hence- 
forth be dropped. The resulting unitary matrix then has 
eigenvalues A± = e*"^ with coscf) = cos0„_q sin {ttv). 

Let us denote N — [W/ L] the integer part of the ratio 
of interface length W and excursion length L = Li = L2- 
Depending on its coordinate y' in the initial Poincare sec- 
tion X — x' , a charge carrier will bounce on the interface 
either TV or iV — 1 times. One can easily see the propor- 
tion of charge carriers bouncing N times on the interface 
is given by the quantity {W/L} = W/L — N. Reflection 
and transmission probabilities for channel {n,a) are then 
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straightforwardly given by the simple expression 

( tz ) - {^/^} ( ) + - ( Zi\ 

(26) 

[This equation is valid actually both in the classical and 
semiclassical frameworks, b ut i n this latter case with = 



Iz^^^P, = \z\^\% From (|25|) we have 



CAT = cos^ {N(j)) + Csin^ 



(27) 




In 



{N(j)) - Csin^ {N(j)) 



with 



sin^ 9n,a — (cos 0„,o, COS (tti/) — sin c 
^ sin^ 9n „ + (cos 0„_o, cos {-kv) — sin < 



(28) 



(29) 



Note that, for a fixed value of the indices {n,a), these 
quantities depend solely on the filling factor u since, by 
the way, so does the angle 9n,a via Eq. (12). 

The total reflection and transmission probabilities R = 
Sn,a ^n,a and T = Y.n,a evaluated in this way as a 
function of interface length are plotted on Fig. |8] It shows 
large oscillations with no sign of emergence of an asymp- 
totic behavior (and of course no conductance plateaus). 
More unexpectedly, the mean value of the semiclassical 
curves in Fig. [8] differs from the classical limit. This is 
also directly visible when comparing classical and semi- 
classical behaviors of individual channels as in Fig. |9] 
It unambiguously signals that interferences between tra- 
jectories at the potential interface dominate the physics 
here. 

The semiclassical behavior can be understood rather 
straightforwardly from the scattering matrix picture. 
The matrix ( 25 1 can indeed be interpreted as that of 



a rotation operator on the Bloch sphere acting on vec- 
tor Zi, defining in this way a discrete map on the 
Bloch sphere. Writing this rotation operator in the form 
-Rn(w) — cos (a;/2)l — isin (a;/2)n • ci, where n is the ro- 
tation axis on the Bloch sphere and ui the rotation angle, 
and identifying with the unitary matrix in Eq. (25), one 
gets u) = 2(1) and 



1 



n 



sm ( 



- sin 9n^a sin (tti^) 
' sin 9n,a cos (ttz^) 

cos 9n^a COS {tTu) 



(30) 



For an incoming electron (whose initial Bloch vector Zq 
points at the north pole), each excursion along the in- 
terface thus amounts to a rotation, on the Bloch sphere, 
of angle w and around the axis given by Eq. ( 30 1 . Two 



limiting cases furthermore provide us with a particularly 
simple picture. Indeed, if the filling factor is an integer, 
i.e. V = Umax [63] we have 





— sin 9n. 

cos 9n,a 



[integer v] 



(31) 



50 

W/1 (interface length) 




W/lg (intetface length) 
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FIG. 8: (Color online) Reflection probability R (top) and 
transmission probability T (bottom) as a function of the 
length of the interface W (in units of the magnetic length 
Ib)- The electrostatic potentials are V2 = — Vi = 3.57 (in 
units of e_B = Hvp/Ib), which corresponds to filling ratios 
1^1 = U2 = (V2/eB)^/2 = 6.37, and thus 13 channels (7 states). 
The mean value of these functions (dotted, black online) dif- 
fers from the classical limit (dashed, red online). The reflec- 
tion and transmission probabilities of each channel behave 
in a similar fashion, oscillating strongly as a function of W. 
The variance of these oscillations is not expected to diminish 
asymptotically as IK — > +00 (see text). 



In that case the axis of rotation depends on 9n,a, be- 
ing close to z for the smallest (0, 1, • • • ) and the largest 
(j^max, "-max— 1, ' ' ' ) channel numbers, and near the equa- 
tor for n ~ 'T.inax/2. On the other hand, the angle of 
rotation is it for everybody, implying in particular that 
there is total reflection for an even number of bounces. 
This is illustrated on Fig. [Toj where we observe that for 
the integer-!/ case considered here, the mean value of the 
transmission and reflection differ significantly from the 
classical 1/2 value. 

If, on the other hand, v = r^max + 1/2 lies midway 
between two integers, we have 



n : 




29n^a [half-integer ly] . (32) 



The axis of rotation is then n 



-X, and is thus within 
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FIG. 9: (Color online) Semiclassical (solid, black online) and 
classical (dasiied, red online) reflection probabilities, as a 
function of W (in units of Ib), for the channel (n, a) — (2, +) 
(which has Klein transmission probability T{62,+) = 0.87). 
The characteristics of the ribbon are the same as in Fig. [8] 
Note that the semiclassical curves are locally straight lines due 
to the piecewise linear form of the reflection and transmission 
probabilities in Eq. (26 1. 



FIG. 11: (Color online) Semiclassical (solid, black online) and 
classical (dashed, red online) reflection probabilities of the 
lowest channel for a half-integer filling factor v — 6.5. This 
case is characterized by an axis of rotation n of iin(i^) lying 
in the equatorial plane and therefore the mean value of the 
reflection probability coincides with the classical limit (to be 
compared with the upper chart of Fig. |10[ ). Also note the 
estimated wavelength of the oscillation between reflection and 
transmission — 'jt\/2u ~ 11 (in units of Ib) is observable 
on this plot. 
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FIG. 10; (Color online) Semiclassical (solid, black online) and 
classical (dashed, red online) reflection probabilities for the 
individual channels (n, a) = (0,—), (2,+) and (5,—). The 
ribbon considered has a filling factor u = 6. The filling factor 
1/ being an integer, this configuration is characterized by total 
reflection of each channel every even bounce and by mean 
values significantly different from the classical limit, especially 
for small or large angles 9n,c,, as is illustrated in top and 
bottom charts. For intermediate angles (middle chart), the 
deviation is less pronounced but still noticeable (compare this 
chart with Fig. ^ . 



the equator and independent on O^.a- As a consequence, 
the mean value of of the transmission and reflection co- 
efficients will correspond in that case (and in that case 
only) to the classical value 1/2. The angle of rotation 
is now however (6'„_o,)-dependent and [modulo {2tt)] is 
close to zero for the smallest and the largest channel 



numbers, and close to tt for n ~ 7i,nax/2. One should 
bear in mind however that Li — 2Ri sin(0„ q,). The total 
rotation ujtot = Nuj, where N ~ W/Li is the number of 
excursions necessary to cross the junction, is thus such 
that 



Wtot 



1 w en,c. 

Phj Ib Sm{dn^a) 



(33) 



For small or large channel numbers, the last factor is 
of order one and thus, if the width of the junction is 
measured in units of Ib, the wavelength of the oscilla- 
tion between transmission and reflection as a function of 
W/Ib is Ai/ — Tr\/2iy, which is indeed what is observed on 
Fig. [IT] This wavelength is reduced by a factor ~ 7r/2 
for intermediate values n ~ iy/2. 

Finally Fig. [12] illustrates an intermediate situation 
{f = riinax + 1/4) for which both n and uj are (0„.a)- 
dependent. We find a semiclassical reflection which is 
in this case lower than in the integer case, but still no- 
ticeably larger than the classical 1/2 value, as well as a 
wavelength for the oscillation between consecutive reflec- 
tions which scales as 27r-\/2z/. 

In summary, the classical transmission through a sym- 
metric n-p junction in the electric regime /3 <C 1 coincides 
with that of Eq. ([T]) . Coherent transport through the n-p 
interface gives rise to interference effects that are charac- 
terized by quantum fluctuations which depend on W/Ib 
and i^. Those fluctuations are different than the UCF 
predicted by RMT. 

The symmetric case corresponds to the special situa- 
tion where the doping is the same at both n and p regions. 
There is no a priori reason to expect the transmission to 
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FIG. 12: (Color online) Semiclassical (solid, black online) and 
classical (dashed, red online) reflection probabilities of the 
lowest channel for a "generic" filling factor (i.e. neither integer 



nor half-integer) i/ = 6.25. Comparing this plot with Fig. 11 
and with the upper chart of Fig. |10| we see that this case is 
an intermediate situation in terms of deviation of the mean 
value from the classical limit. 



be the same for both V2 = —V\ and V2 7^ —V\. Hence, to 
make contact with experiments, we need a theory for the 
general case where the gate voltages are not symmetric. 
This is what we do next. 



V. ELECTRIC REGIME: GENERAL CASE 

This section is devoted to the calculation of the trans- 
mission through a n-p junction for the general case, where 
V\ 7^ — V2, in the electric regime. For that purpose, it is 
no longer possible to use the intuitive scattering matrix 
approach discussed in the previous section, which is valid 
for the symmetric case only. In what follows we develop 
a semiclassical theory that allows the calculation of the 
Landauer transmission for the general case. We start 
discussing how to apply the Fisher-Lee formula (SO] for 
a graphene n--p junction. Then we address the dynamics 
at the interface, which provides the elements required by 
the semiclassical calculation. We proceed presenting the 
main technical details of the derivation, complemented 
by some additional material presented in the appendices. 
We conclude this section discussing the classical limit and 
summarizing the main features of the transmission in the 
general case. 



A. Fisher-Lee/Baranger-Stone formalism 

We now consider the electric regime of a step- like junc- 
tion for arbitrary values of V\ and V2 (V2V1 < 0), for 
which we obtain a semiclassical evaluation of the Lan- 
dauer conductance, Eq. ([2|. This is achieved with the 
help of a formalism which was first introduced by Fisher 



and Lee [50] j and later generalized to account for a mag- 
netic field by Baranger and Stone [3T]. This formalism 
is based on the use of Green's functions for which we 
derived in a previous work S5] a semiclassical approxi- 
mation in graphene. As was discussed in that paper, the 
distinguishing feature of this Green's function as com- 
pared with the standard 2DEG Schrodinger expression is 
the appearance of a semiclassical phase which can be un- 
derstood as the topological part of the usual Berry phase 
occurring in the context of systems depending adiabati- 
cally on an external parameter. However in the absence 
of a mass term in the graphene Hamiltonian (as will be 
the case in this work), both phases are equal and can 
be expressed as i^c = - 6'')/2, with B' and Q" the 

angle of the initial and final pseudo-momentum H of the 
corresponding trajectory. We shall thus refer to it in the 
following as the Berry phase. 

Turning back to the Fisher-Lee/Baranger-Stone for- 
malism, the formulae obtained in [Ml EI] were derived 
for Schrodinger (scalar) electrons and should be some- 
what modified to describe charge carriers in graphene. 
Special attention must be paid to the pseudo-spin degree 
of freedom which shows up in the spinor structure of the 
modes Xn,a and the matrix structure of the Green's func- 
tion and which generates non-commutative operations. 
With this in mind, calculations are rather straightfor- 
ward and the following expressions can be obtained for 
a general mesoscopic graphene sample with an arbitrary 
number of leads (/i,J^,- • • ): the conductance from lead v 
to lead \i (with /i 7^ ;/) is 

o2f 



(34) 



X ^(r,„r',)(VpH-e,)g^(r^,r', 



and the transmission probability amplitude of going from 
channel n in lead v to channel m in lead fi is 



f — —ih 



C„ JC. 



dVf, / dyUm (I'm) (Vp'H • e^) 



(35) 



With these notations, the Landauer-Biittiker formula 
reads 5^1, = X]„ Z^m ^nd are trans- 

verse sections of the leads, while and e,y stand for 
the unit normal (outward pointing) vectors to the corre- 
sponding sections. "H is the graphene Hamiltonian, Q the 
(retarded) Green's function and Xn the quantized mode 
(with ± labelling its direction of propagation with respect 
to the central region separating the leads). Dependence 
on the valley index a in the modes has been temporarily 
dropped for convenience. 



The expressions ( 34 ) and ( 35 1 can be slightly lightened 
when the two valleys K and K' are uncoupled and can be 
treated independently. Starting from the valley isotropic 
representation ^ = {iI;a,'4'b,—'4''bt'4''a)'^ introduced in 



section |HH the effective Hamiltonian within the valley a 
can be written 



H = vpU ■ a + U{r)l + am{r)vpaz 



(36) 
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with the convention that 



* = 



i>2 



i-i^'s^^'^V if a = -l 



(37) 



With this choice of representation, expressions (34) and 



(35) read 



li^t" = -G, 



/ dy^, / dy'^Tr[{a ■ e^) 
Jc. Jc„ 



-ihvp 



dy, 



(38) 



(39) 



Focusing now on the specific geometry under consider- 
ation (cf. Fig. [6]), we can drop the lead indices, assume 
the sections from which the conductance is computed to 
be located at the extremities of the junction (at abscissa 
x' on the incoming side and x" = x' + on the outgoing 
one) and use the coordinate y inside the section. The 
transmission coefficients can then be written as 

trnn^ihvl [ dy" [ dy' x+H^'M {v" , r')a,Xn i^') ■ 



(40) 

Note that since we are working in the quantum Hall 
regime W ^ Ib, integrals in Eq. (40) are effectively re- 
stricted to one of the edges. 

As we are interested in the total conductance rather 
than the individual transmission coefficients, we do not 
need to compute all the t-mn but only the sum J2m KmnP- 
Using that 



dyx^\r)vFcr^Xni^) = 



(41) 



as is proven in 31j, one can easily show that 

T„,„ "i/ J2 l^""!' = / dy"TL{r")vFa,% 
Jv">o 

with 

r„,a(r") =ih f dy'g{r",r')vF<J.XnJ.r') ■ (43) 



(42) 



The same expressions apply for Rn.a = 1 ^ Tn.a, except 
that in Eq. ( 42 ) the integral should be taken in the elec- 



tron side of the junction, i.e. on y" < 0. The prescription 
of directly computing T„^ct instead of the individual tmn 
additionally bypasses the need to project the incoming 
modes Xn propagated along the interface on the outgo- 
ing ones Xm iu Eq. (40). This is particularly useful on 



the p side of the junction where angle 02 has no reason 
to coincide with a quantized value and where transmit- 
ted charge carriers are therefore no longer in a properly 
quantized state but in a superposition of outgoing modes. 



Our main task is now to evaluate semiclassically 
Eq. ( 43 ) . This requires obtaining semiclassical approxi- 
mations of the incoming mode X^c(r') and of the Green's 
function Q{r",r'). The mode Xna(i"') is built semiclas- 
sically on the manifold obtained from the one-parameter 
family of trajectories bouncing with an angle 9n,a on the 
edge of the lead. Within the representation (37) and 



sticking with an infinite mass edge confinement, one gets 



Xn,Q (^) 



\/|sin6'„,a(y)| ^f^^ 



g-j§e,.,a(a) 



e'2 



(44) 



with = ± an index (not to be confused with the filling 
factors i^i) labeling the sheets of the manifold on which 
the mode is constructed ( Pj, > for ;/ = -1-1, pj, < for 
v — —1). The caustic in phase space at the junction of 
the two sheets is taken into account by the phase jump 
/i(j^) = 0(— i^) (with 9 the Heaviside step function). In 
Eq. (44) 9n,a{y) = cos^^ {cos9n,a — V / Ri) is the angle of 



the tangent to the trajectory when at a distance y from 
the edge, A:"'" = kp cos 9n,a is the constant of motion 
associated with the mode, Sn,a{y) = ^^ifi9n,aiy)) the 
action, and Cn,a = (4t'_F-Ri sin0„_a)~^/^ is a normaliza- 



tion constant which is determined from Eq. (41 ) 



Turning now to the Green's function, a semiclassical 
approximation valid in either the electron or hole re- 
gion was derived in [32]. Including Klein tunneling, i.e. 
the transitions from electron to hole regions to this for- 
malism, is a priori a non-trivial (although feasible) task 
in a completely general setup. The limit /3 ^ 1 that 
we consider here, and the fact that we assume a per- 
fectly straight potential step, simplify however consider- 
ably the problem. Indeed one can in this case, for the 
Klein tunneling, treat the semiclassical wavefunctions as 
plane waves, and therefore use the transmission and re- 
flection coefficients Eqs. (18)-(19) as in [JB]. This leads 



to the straightforward generalization of the semiclassical 
Green's function expression 

g..(r",r';.,) = ^[nn^^"^^i''') 

j \i=l-fi,rii / ^^g^ 



in^/2TTin\Jj{r",r')\ 



Vf''"\r")vf'''^\r') , 



with 7.i, rji labelling scattering events on the n-p junction, 
respectively associated to reflections and transmissions 
in region i. As usual, the sum is over all classical tra- 
jectories joining points r' and r" at the Fermi energy ep. 
The phases accumulated along the way include the action 
Sj , the Maslov index /i^ counting the number of caustics 
traversed, and the Berry phase — ^{9" — 9j)/2, with 
9j and 9'J the direction of the initial and final pseudo- 
momentum n of the trajectory j. Note that as the total 
Berry phase is accounted for in , the reduced refiection 
and transmission coefficients (without the Berry phases) 
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ri and ti given by (18)-(19) should be used in Eq. (45 1 
The determinant 



n = 3 



J, 



dy"dy' 



(46) 



implements the conservation of classical probability. Fi- 
nally, V^'''^\r) is the eigenstate along the j*'* trajectory 

of the classical Hamiltonian i?''^''' = V{x) + e{r)vp\'n\, 
with e(r) = ±1 (+1 on the electron side and —1 on 
the hole side). In the absence of any mass term in 
the bulk of the sample and choosing the representation 
(37), these eigenstates depend solely on the angle of 
the pseudo-momentum: Vj^{r) = (l/\/2)(l, e*^J )"^ and 

Determining the specific Green's function for the prob- 
lem under consideration can essentially be reduced to the 
task of making a complete list of the trajectories connect- 
ing the Poincare sections on both sides of the interface 
and computing the probability amplitudes associated to 
each one of them. This issue will now be addressed. 



B. Dynamics at the interface 

Depending on the relative size of the potentials |Vi| 
and V21 Eq. (15) defines a critical angle for either or 



Bi above which reflection on the n-p junction is total. 
Without loss of generality, we will take < V2 which 
constrains incident holes on the interface to the angular 
domain [6'crit,7r-6'c„t], with 6*^^* = cos"^ (|V^i|/V^2)- This 
allows not to worry about possible total reflection on the 
electron side, and additionally sets the excursion length 
scales L2 > Li. 

Consider a typical trajectory going from r' = {x',y') 
in the initial x — x' Poincare section to r" = {x",y") in 
the final x = x" = x' + W Poincare section. The trajec- 
tory can be labeled by an index j specifying whether the 
trajectory is transmitted or reflected for each of its suc- 
cessive encounters with the junction. For a given j, the 
initial and final angles 6' and 9" are fixed once the coor- 
dinates y' and y" are. The trajectory j, that is the suc- 
cessive list of transmissions and reflections, can be char- 
acterized by two integers: the number of excursions n2 in 
region 2 (which then fixes the remaining number ^1(712) 
of excursions in region 1) and the number of traversals k 
of the interface. A typical example is shown in Fig. [13] 

These two integers do not specify uniquely the trajec- 
tory j since it is possible to permute the order of the ex- 
cursions in the hole and electron regions while preserving 
the couple {n2,k). They however define a class of trajec- 
tories which, as we will see, give the same contribution to 
the semiclassical Green's function Eq. (45). Indeed, the 
mapping {y',0') {y",9") depends only on n2, which 
fixes the number of excursions in regions 1 and 2, and on 
the parity of k, which determines whether the trajectory 
j exits the junction in region 1 (fc even, reflection) or 2 (fc 




FIG. 13: (Color online) Trajectory characterized by n2 = 
3 excursions on the hole side and k — 4 traversals of the 
interface 



odd, transmission). This mapping remains however un- 
changed if the order of the excursions in regions 1 and 2 
is modified. As a consequence the determinant Eq. ( 46 1 , 



which can be expressed in terms of this mapping, or the 
Berry phase = —{9" — 9')/2 which is only a function 
of the initial and final angles of the trajectory, are also 
independent of the ordering of the excursions. 

In the same way, the action Sj and the Maslov in- 
dex f i j i n the phase of the semiclassical Green's function 
Eq. ( 45 ) are functions of n2 and the parity of k only. 



The action, for instance, can be expressed as Sj{n2) = 
hk^W + ni{n2)Si+n2S2 + /^S, with Si = 2hvif{9i) and 
^2 — —2hv2f{'K — O2)) the actions accumulated along 
an excursion in regions 1 and 2 respectively (note that 
52 < since, in region 2, p • r < 0), and 

AS = hui {f{9i)-f{9')) 

hui {f{9") - f{-9i)) [k even, reflection) 

fii^2 {f{0") — /(tt — 92)) (k odd, transmission) 

(47) 



+ 



9^9— (sin 26*) /2 is the same function as in sec- 
III). The angles 9i and 6*2 are the ones introduced 



(/ : 
tion 

in Fig. |6J and at this point should be understood as be- 
ing functions of 9' and 9" . Explicit computation of the 
Maslov index Hj{n2) (see appendix IaI) shows also that, 
quite naturally, it does not depend either on the ordering 
of the excursions. 

We now turn our attention to the factors associated 
with scattering at the interface. Let us first consider the 
case where the trajectory exits the junction in region 1, 
i.e. of an even number of traversals k = 2k'. In that 
case, the number of traversals from 1 to 2 as well as from 
2 to 1 are equal to k' , and there are n2 — k' reflections 
on side 2 and ni(n2) -|- 1 — fc' reflections on side 1 (the 
additional term 1 coming from the fact that the trajec- 
tory initially leaves the Poincare section in that region). 
The probability amplitude associated with the reflections 
and transmissions at the junction interface for the class 
of trajectories (712, fc') is thus given by 



AR{n2,k') 



«l(n2) + l n2 
'1 '2 



^1^2 
7-1^2 



(48) 



The case of an odd number of traversals k = 2k" + 1, 
i.e. when the trajectory j is eventually transmitted in 
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FIG. 14: (Color online) Class of trajectories characterized by 
k' = 3 (and an arbitrary n2). Filled (blue online) rectangles 
symbolize a succession of excursions in region 1 and hatched 
(red online) rectangles a succession of excursions in region 2. 
Although these rectangles are represented as if they had the 
same size, each one of them may contain a different number 
of excursions. 



region 2, is completely equivalent. There are then k" + 1 
transmissions from region 1 to region 2, k" transmissions 
from region 2 to region 1, n2 — k" reflections in region 2 
and ni(n2) — fc" in region 1. The probability amplitude 
for the class (n2,fc") thus reads 

^T(n2,fc") = tirr^"=^V2"^ (^^y . (49) 

All permutations in the order of the excursions along 
the interface that preserve numbers n2 and k correspond 
to the same probability amplitude. We must therefore 
determine the degeneracy factor f2 that gives the number 
of distinct trajectories belonging to the class character- 
ized by integers (n2,fc). Starting again with k = 2k' even 
(reflection) , let us materialize each uninterrupted succes- 
sion of excursions in region 1 or 2 by a rectangle (see 
Fig. |14[). If we include excursion portions respectively 



leaving from the initial Poincare section and arriving at 
the final Poincare section, there are fc' rectangles in re- 
gion 2 and k' + 1 rectangles in region 1, when the total 
number of excursions in region 2 is rt2, and that in region 
1 is 711(712) -l-2. Then, using the combinatorial result that 

there are ( ^ ways of writing an integer 77 as a sum 

of k non-zero integers (or equivalently of distributing n 
excursions into k non-empty intervals), the degeneracy 
factor is straightforwardly given by 

(SO) 

with = k\{n-ky. binomial coefficient. 

For an odd number k — 2k" -|- 1 of traversals (trans- 
mission), there should now be an equal number k" -\- 1 
of rectangles in regions 1 and 2, while the total number 
of excursions are respectively 711(772) -I- 1 in region 1 and 
772 + 1 in region 2, once more including here "partial" 
excursions leaving from the initial Poincare section and 
arriving at the final Poincare section. Using the same 
combinatorial result as before, this yields for the degen- 
eracy factor 

o.,-..r,.(-)(-<-))^ (-) 

Combining all of the results obtained in this subsec- 
tion, the semiclassical Green's function for our problem 
can be expressed as 



g«,(r",r')= y : ^ ^ ^^ .^(„,).tg("2)-^fM(».)+.C 
f-„7V2^7;7|J„,(r",r')| 



(52) 



with 



S(772) 



Y.^, ^fl(772, k')^R{n2, k') if y" < (refiection) 
^j,„ At(772, k")nT{n2,k") if y" > (transmission) 



(53) 



and N2 the upper bound of the number of excursions in 
region 2 (bounds for the number of traversals are given 
by min(l,772) < fc' < min(772, 771(77.2) -|- 1) and < A:" < 
min(772, 771(772)). 



C. Semiclassical expression for the conductance 

Let us first discuss the case of reflected trajectories 
{y" < 0, e(r") > 0). Using the semiclassical expressions 



( 44 1 and ( 52 ) , the matrix structure in integral ( 43 ) reads 



X COS 



(54) 



Computing the integral ( |43[ ) in the semiclassi- 
cal limit S — > can be done using the sta- 
tionary phase approximation J A{x)ei^^^'^dx ~ 
^(x,)eK^'(^=)(27ri;7/|525'(x,)|)i/2e-*fA's^ with 
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Xs the points where S{x) is stationary and 
^jLs = &[—d'^S/dx^{xs)]- The stationary phase con- 
dition 



dy' 



= 



^'Or.Ay's) = o'M) (55) 



V",n2 



expresses that for a given final position y" (and a given 
number of excursions in region 2), the stationary phase 
point j/j is the one where the initial angle matches with 



that of the quantized mode. This implies the additional 
identifications 9i = dn,a, H = k^'" and 9'J = v"9n,a{y") 
with v" the sheet index in the final Poincare section. The 
action of the Green's function thus becomes Sj^y'^) = 

hk^^'^ix" - X') + (ni + l)5f'" + 712^2 + v"Sr,,o.{y") - 

v'Sna{y's)j with 5"'" ~ 2hvif{6n.a) quantized as in 
Eq. |8]). Inserting these results in the integral (43), and 



bearing in mind that the stationary phase point y'g de- 
pends on the integer n2 we obtain 



v"=±l n2=0 



dl,{S, + SnA (2/;)J,(2/;)sin0„,„(yO 



fe' 



(56) 



In Eq. (56), ^ij^s is the sum of the Maslov index 
in the Green's function and of the index jj,s coming 
from the stationary phase integral. The latter is zero if 
dl,{Sj + Sn,a) > and one if dl,{Sj + Sn^ < 0, while 
the former requires some care to be computed precisely. 
The technical calculation of /ij is detailed in appendix 
[A] For our current purposes, one can actually show that 
^J'J,S = K^') - K^") + ni + 1 - n2 with ^i{v") = 
the phase jump at the caustic between the sheets in the 
final Poincare section. 

The final step of this calculation involves computing 
the prefactor (9^,(5',- -I- S'„,q)| „ J,)^^/^, which we do in 



appendix B Inserting the result in Eq. (56), one finds 



that all trace of the stationary phase point y'^ has van- 
ished and that expression ( 56 ) can be simply written as 



N2 



n2=0 
X 



(57) 



YAR{n2,k')^R{n2,k') 



Comparing this with the original integral (43) makes it 



possible to give a rather transparent interpretation for 
the role played by the Green's function. It basically 
amounts to propagating the original mode from r' to r" 
with a certain probability weight corresponding to the 
various trajectories fulfilling the stationary phase con- 
dition (55 1 and connecting these points. The reflection 



probability for channel n polarized in valley a is obtained 
by inserting Eq. ( 57 ) in Eq. ( 42 ) , which immediately gives 



Rum = VF\Cn,a\^ J dy" 



2 COS 6" 



smf 



JV2 



J2 (-zei«"'°)"^+i(zei^^)"^ ^ Anin^, fc')^^/?(^2, k') 

n,2— k' 

(58) 



with the shorthand 6" — 0n,a{y")- A change of variables 
dy" — Ri sin 9"d6" leads to the final result 



-ST 



Rn,a = / dO' 



cos ( 



9^ 2sin6'„,Q 



N2{9") 



Y (-^e^^"'°)"l("^■''")+l(^e^^^)"^^AflX^2,fc')^^i?.(^2,fc') 



ri2=0 



(59) 



with the integer-valued functions N2 and ni given by 
'W- Ri sm0n.a - Ri sin 9" 



N2{e") 



2i?2 sin 02 



ni{n2,0") 



(60) 



W - Ri sin0n,a~Ri sine" 
2i?i sin On^a 

R2 sin 6*2 



-n2 



i?i sin 0n 



(61) 
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FIG. 15: (Color online) Portion of cyclotronic orbit illustrat- 
ing that, when 6n,a > f , accessible values for 6" do not any 



more correspond to the range [~6 



Note the bound in the integral (59l is 6'^ = mm{9n,a, t^ — 
(^n,a) rather than simply 9n,a- TEs is because as the edge 
angle 6n,a exceeds 7r/2, trajectories with a final angle 9" 
larger than tt — 9n^a traverse the Poincare section twice 
and scatter once more on the interface (see Fig. 15). The 
limiting angle in Eq. ( 59 1 is then the one for which no 



further scattering on the potential step can take place. 
The integral in Eq. ( 59 1 can be estimated numerically as 



a function of the interface length W and the tunable field 
strengths Vi, V2 and B. 

The transmission T„.q can be calculated following sim- 
ilar steps. For the reader interested in the technical de- 
tails, a summary of the derivation is presented in Ap- 
pendix [C] The result reads 



T — 



d9'- 



R2 cos 9" 
1R\ sin 0„ Q 



N2(e") 

E 

"2=0 



(-zes^?r'°)"i("2^»")(ie^'^^)"^ y ATin2,k")nT{n2,k") 



k" 



(62) 



with 9t — min(02, vr — 02)- The integer valued functions 
N2 and ni are given by the formulae 



D. Comparison of classical/semiclassical 
predictions and summary 



N2{9") 



ni{n2,9") 



W - R2 sin( 



R2 sine" 



2i?2 sin 6*2 



'W-R2 sm92 + R2 sin 6*" 



(63) 



2Ri sin 



i?2sin6'2 

—n2 : — 

i?i sin 9n^a 



(64) 



Bounds in the integral (62 1 depend on the sign of 6*2 — 7r/2 



for the same reason as for reflected trajectories. 



Equations ( 59 ) and ( 62 1 readily give the total reflection 
and transmission coefhcients, namely, R = ^„ ^ Rn,a 
and T = J2n a'^n,a- Alternatively, these results can be 
inserted into the Landauer formula, Eq. ([2]), giving the 
conductance. These are, from the technical point of view, 
the main results of this paper. 

In Fig. [16) wc show R and T for a couple of values of 
the electrostatic potentials Vi and V2 ■ The reflection and 
transmission coefficients show an oscillating behavior as 
a function of the interface length W. 

The overall behavior of R and T is qualitatively similar, 
but not identical, to the one observed in the symmetric 
case studied in the previous section. In particular, no 
saturation of the conductance is observed or predicted. 
Before discussing these results with greater depth, let us 
first gain some insight on what is expected classically in 
the general situation \Vi \ 7^ V2. 



Contrary to the semiclassical probabilities derived in 
the previous subsection, their classical counterparts con- 
verge to an asymptotic value in the limit of a long enough 
interface. This can be easily shown using the following 
line of reasoning. Let us call Vi (x) and 7^2 {x) the classi- 
cal probabilities for a charge carrier to be found in region 
1 and region 2 at a longitudinal distance x from the initial 
Poincare section. These probabilities obey equations 



Vi{x) = (1 - T)Viix - Li) + TV2{x - L2) 
V2{x) = (1 - T)V2{x - L2) + TVi{x - Li) 



(65) 



whose only asymptotically constant solution as a; 3> L1.2 
is 7^1 = 7^2 = 1/2. Noticing additionally that charge 
carriers emerging in region 1 (respectively region 2) must 
do so at a distance smaller than Li (respectively smaller 
than L2) from the final Poincare section, one gets, in the 
asymptotic limit, the classical reflection and transmission 
probabilities 



R-l^^L,/{L^+L2) 
T^'^L2/{L^+L2) 



as W 



+00 . 



(66) 



These values are indeed those observed if one plots the 
classical counterparts of the Fisher- Lee formulae ( 59 ) and 



(62 1, which can be obtained by considering electrons and 
holes as non-interfering classical particles. This essen- 
tially amounts to replacing probability amplitudes by 
probabilities and neglecting all phase factors accumu- 
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FIG. 16: (Color online) Total reflection and transmission 
probabilities as a function of W/Ib- Top: for filling fac- 
tors = 6.37 and ^2 ~ 12.15 in the electron and hole re- 
gions respectively. Bottom: for filling factors vi = 3.70 and 
V2 = 6.99. The qualitative behavior of individual channel 
reflection and transmission probabilities is very similar. 



lated along the trajectories, giving 



ncl _ 



COS 0" 
2 sin 9n.a 



N2{9") 



(n2,0")+"2 + l 



n2=0 



T 



T 



2k' 



^R{n2, k') , 



O'j' 



d6 



„ i?2 COS ( 



N2{e") 



2RiS\a.9n,a 



(67) 



)+n2 



112=0 



T 



T 



2k" 



nT{n2,k") 
(68) 



with T = Tie^.o.) = sin0„,„sin(?2/cos2((0„^„-02)/2) 
the Klein transmission probability. These classical for- 
mulae are plotted numerically as a function of the in- 



FIG. 17: (Color online) Classical reflection probabilities of 
individual channels for a filling factor u\ = 6.37 in the 
electron region and V2 ~ 12.15 in the hole region. Each 
channel is characterized by a Klein transmission probability 

T{er,,c): T{eo,-) = 0.45, T{e2,+) = 0.90, r(e»3,+) = 0.991, 

T{9^^-) = 0.78. When T ~ 1, classical convergence is masked 
by large oscillations on the scale of the considered interface 
length W. 



terface length W on Fig. 17 As can be seen, some of 



the probabilities converge rather quickly to the asymp- 
totic values mentioned above, while others show slow con- 
vergence, sometimes barely visible on the scale (value 
of W) used. This is of course simply due to the fact 
that the convergence speed depends on the value of the 
Klein transmission probability. When the latter ap- 
proaches unity, the potential barrier becomes transparent 
for charge carriers which are thus alternatively reflected 
or transmitted. 

Quasi-unit Klein transmission probabilities are ob- 
tained for angles close to 7r/2 which, as illustrated in 
Fig. 18 are more densely sampled when quantizing the 
dispersion relation ( 12 1. This behavior is however mainly 



due to the assumption we made of considering an ex- 
tremely abrupt potential step on the scale of the magnetic 
length. Restoring a finite steepness to the n-p junction 



would have the dashed (green online) curve in Fig. 18 
look much more like a sharp peak and considerably re- 
duce the likeliness of having quantized angles with a close 
to unit Klein transmission probability [64j . 



Coming back to the semiclassical formulae ( 59 1 and 



(62) plotted in Fig. 16 and as we already pointed out. 



their behavior is qualitatively very similar to what was 
observed in the symmetric case plotted in Fig. [S] We 
lack here the equivalent of the Bloch sphere picture valid 
in the symmetric case, but we believe nothing funda- 
mentally different is going on here and the physics is 
hence essentially the same. The presence of large os- 
cillations in the semiclassical transmissions, as opposed 
to the classical ones, once again indicates that interfer- 
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6/71 (angle of incidence) 



FIG. 18: (Color online) Functions f{0)/-K (solid, black online) 
and T{6) (dashed, green online). Squares (red online) mark 
the values of the quantized angles for the filling factor i/i = 
6.37. As can be clearly seen, most quantized angles sit in the 
central region 9 £ [n /3,2n /3], bounded by the vertical dot- 
dashed lines, and three of them are sufficiently near n/2 that 
the semiclassical behavior and the classical one are almost 
indistinguishable for the interface lengths considered. 



ences between trajectories at the potential interface are 
strong. Concerning the mean values of the probabilities 
plotted in Fig. 16 we found that: (i) They still differ 
from the classical limit (as in the symmetric case), and 



(ii) They compare poorly with the full mode mixing hy- 
pothesis prediction, given by Eq. ([T]) . These discrepancies 
clearly invalidate the possibility of full chaotic mixing in 
our model, which can be straightforwardly understood 
by observing that the accessible region of phase space 
on the hole-doped side of the junction is restricted by 
Eq. ( |15[ ). A better approximation to the mean semiclas- 
sical value can actually be obtained using a mode mixing 
hypothesis with an effective number of modes in region 
2 corresponding to the accessible region of phase space 
(angles 6*2 in [-Ocrit^Ocnt]), i-e. 



N2 ^ N!, 



, def 1 



Pydy ~ N2II 



fide 



... 

Partitioning the current with equal probability in the 
available edge channels yields a reasonable approxima- 
tion (about 10%) to the mean conductance: {G/Go)w — 
NiNy{Ni+N^). 

Summarizing, in this section we show that, in the semi- 
classical limit, the general (non-symmetric) case shows 
large transmission/reflection fluctuations. Furthermore, 
suppressing these fluctuations by averaging over W would 
not recover the full mode mixing prediction. A reason- 
ably good approximation of the mean transmission or re- 
flection values can however be obtained from a modified 
mode mixing hypothesis where only the number of modes 
corresponding to the accessible phase space is taken into 



account. 

As in the symmetric case, mesoscopic corrections lead 
unavoidably to quantum fluctuations. Those are clearly 
non-universal: for certain values of Vi and V2 the magni- 
tude of the fluctuations decreases with increasing W/l b , 
while for other combinations of Vi and V2 the fluctua- 
tions do not seem to depend on W/l b and no systematic 
behavior is observed. Notice that we investigate W/Ib 
values that are comparable with the experimental ones. 
The variety of fluctuation patterns we observe can be 
semiclassically explained by means of the quantum inter- 
ference between different snake-like trajectories. 



VI. MODEL ASSUMPTIONS, LIMITATIONS, 
AND EXTENSIONS 



The purpose of this section is to discuss the most im- 
portant assumptions made in our model. Some of them 
were introduced for the sake of simplicity and do not 
introduce limitations to our analysis. Others were nec- 
essary to proceed analytically and need further justiflca- 
tion. 

We begin addressing the geometry of our transport 
conflguration. Recall that instead of the experimental 
setup of a perpendicular junction, we have settled for the 
analytically simpler setup of a "parallel" junction (see the 
discussion at the beginning of section IV). This issue is 
addressed in Ref. [45] where we demonstrate that, aside 
for possible diffraction effects at the edge-junction cor- 
ner, and despite a notable difference in complexity when 
it comes to exhaustively identifying all trajectories con- 
necting the initial and final Poincare sections, the con- 
ductance is very similar in both configurations. 

Let us now discuss the effect of different boundary con- 
ditions on the conductance. For the magnetic regime, 
this has already been discussed at the end of section [Till 
so we shall limit the discussion here to the electric regime. 
The main effect of choosing zigzag or armchair edges in- 
stead of an infinite mass confinement is to modify the 
quantization condition Q and thereby the values of the 
quantized angles for the edge channels. The EBK semi- 
classical quantization procedure was recently applied to 
the zigzag and armchair cases |41) and yields quantiza- 
tion conditions which can be obtained very simply from 
the infinite mass case by computing the new phase shift 
(or (j)a) acquired by a plane wave scattering on a zigzag 
(or armchair) edge in graphene. As we have seen how- 
ever, the conductance of individual channels did not show 
any special feature dependent on the value of the quan- 
tized edge angles 0„_q. All edge channels had a qualita- 
tively similar behavior, their conductance oscillating as 
a function of W. We therefore do not expect that taking 
zigzag or armchair boundary conditions (which basically 
amounts to changing the quantized values of the edge 
angles) will qualitatively alter our predictions. 

We now address the effect of replacing our step-like in- 
terface by one with a finite steepness L <ti Ib- As long 
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as electrons can locally still be approximated as plane 
waves near the interface, a finite steepness in the po- 
tential step essentially sharpens significantly the angular 
profile of the Klein transmission probability (the dashed 
green curve in Fig. 18 ) around the angle of perfect trans- 
mission. An analytical expression for the Klein trans- 
mission probability in this context was derived by Shy- 
tov and collaborators |33j. It is found that the Klein 
transmission probability is perfect for the incident angle 
Op = TT — cos~^/3 [65, and decreases exponentially with 
L as the angle of incidence is brought away from 9 p. A 
slight asymmetry between Klein transmission probabil- 
ities on both sides of the junction is additionally cre- 
ated when —Vi ^ ¥2- Nevertheless, the results presented 
in this paper should qualitatively hold true as long as 
L <^ Ib (or equivalently /3 ^ 1). 

Let us end our "tour d'horizon" by discussing the ef- 
fect of disorder in the system. Decreasing further the 
steepness of the barrier favors the charge carriers to dwell 
longer at the vicinity of the Dirac point which, in most 
experimental setups, is characterized by the existence of 
electron-hole puddles [17] combined with a weak screen- 
ing of electrostatic charges (due to the vanishing den- 
sity of states). These effects tend to enhance the influ- 
ence of impurities on the electrons at the interface region 
PSI , and possibly drive the system out of the fully 
coherent ballistic regime that we consider in this paper. 
Evidently, the inclusion of disorder in our model favors 
a transition to the regime described by RMT. Indeed, 
the presence of impurities at the junction interface would 
randomize the scattering angles and suppress the restric- 
tion on the available phase space in region 2 imposed by 
Eq. (15 1, bringing the average conductance closer to the 
full mode mixing value, Eq. Q. We however expect a 
rather smooth transition from the ballistic non-universal 
fluctuations we calculate here to the universal ones, ubiq- 
uitous in chaotic and disordered systems. 



VII. CONCLUSION 

We have studied electronic transport in a graphene n- 
p junction subjected to a strong perpendicular magnetic 
field. Our main interest in this problem was twofold. 
Our first goal was to shed some light on the experimen- 
tally observed conductance plateaus in this configuration, 
which still lack a full theoretical explanation. Second, 
we wanted to confront the full mode mixing hypothesis 
with a full analysis of a model as consistent as possible 
with the most physically relevant parameters. The lat- 
ter was further motivated by the fact that an elementary 
classical calculation predicts equal reflection and trans- 
mission probabilities through a symmetric graphene n-p 
junction in the Hall regime, which raised the question 
of how much this property would remain true within a 
semiclassical description. 

Concerning our first goal, in distinction to the UCF ob- 
served in numerical simulations of disordered n-p junc- 



tions [251 J we obtain large non-universal transmis- 
sion fluctuations in ballistic junctions. The combina- 
tion of these results essentially rules out the possibility 
of explaining the observations of fractional conductance 
plateaus within a fully coherent description of the junc- 
tion. To reconcile experimental and theoretical results, 
we believe that a decoherence mechanism suppressing the 
interference effects is necessary |20| . Among the various 
possibilities, such as interactions with electron-hole pud- 
dles, electron-phonon and electron-electron interactions, 
one should find one that is particularly effective in the 
vicinity of the junction. For instance, it is plausible that 
experimental n-p junctions are not as abrupt as the ones 
considered in our model, belonging to an intermediate 
regime of ~ 1 where a random network of electron- 
hole puddles |50j could provide both random current 
partitioning and decoherence. With the recent advent 
of hBN substrates [5TH53] which, when intercalated be- 
tween graphene and Si02, were shown to significantly 
increase electrical mobilities and equally reduce charge 
density inhomogeneities, we hope new experiments can 
shed more light on the nature of the dephasing mecha- 
nism taking place at the junction interface. The effect of 
inelastic scattering events near the Dirac point certainly 
also deserves future investigation [54]. 

Concerning our second goal, as already stated, we ob- 
serve that interference effects play a dramatic role for the 
transport in ballistic n-p junctions irrespective of the ra- 
tio W/l B ■ Our expectation of a self-averaging mechanism 
in the summation over a large number of semiclassical 
contributing terms was not fulfilled. Surprisingly, in the 
general case of Vl 7^ — V2, the mean transport quanti- 
ties (reflection and transmission probabilities) expected 
from the classical dynamics of the junction differ consid- 
erably from their semiclassical (and thus quantum) coun- 
terparts. In the case of symmetric junctions, for which a 
simple and transparent scattering matrix approach can 
be used, these interference effects can furthermore be de- 
scribed within a Bloch sphere picture, which makes it 
possible to give a natural interpretation of the discrep- 
ancies between classical and semiclassical behaviors. 

Our results suggest the need of further experimental 
insight to understand the absence of transmission fluctu- 
ations in n-p junctions at the quantum Hall regime. We 
believe that the conductance fluctuations predicted here 
should be observable in junctions for which the decoher- 
ence processes are suppressed, particularly for suspended 
graphene. 
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FIG. 19: (Color online) Top: Caustics along a reflected tra- 
jectory with (Ji < 0. The rightmost caustic concerns the final 
portion L" and originates from the fact that 02 > 0. Bot- 
tom: Variation of the number of caustics /ij as ai becomes 
positive. If ^1 < n/2, the transitional excursion contains no 
caustic, while if 61 > ■k/2 the transitional excursion contains 
two caustics. 



Appendix A: Calculation of the Maslov index /xj 



In this appendix we compute the Maslov indices that 



appear in Eqs. (57) and (C2). Recall the Green's func- 



tion Maslov index fij counts the number of caustics on 
the trajectory j, a caustic being a point where j and 
the neighboring trajectories obtained by an infinitesi- 
mal change of the initial momentum intersect each other. 
Here however it should be noted that, since in the hole re- 
gion momentum and velocity have an opposite direction, 
the caustics there should be counted negatively. 

Caustics can be found in three different parts of the 
trajectory: before the first encounter of the trajectory 
with the interface of the junction; after the last encounter 
with the interface; and in-between these two points. We 
shall call respectively /ij, /i" and /ij the three corre- 
sponding contributions to the Maslov index, with /ij = 
fj,j + ^ij + fj,j. We further note L' = i?i(sin0i — sin6'') 
the distance between the first encounter and the initial 
Poincare section, and L" = Ri{s'm6" + sin^i) {i = 1 
for reflected trajectories, i — 2 for transmitted trajecto- 
ries) the distance between the last encounter and the final 
Poincare section. Reflected trajectories and transmitted 
ones will be addressed separately. 



1. Reflected trajectories 



Let us introduce the following indices 



(Ti = sign 



(72 = sign 
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sm 
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_i dL" 



cos 6*' - (2ni -f 1) cos 6*1- 



smf 



v' 



-2ri2 cos 92- 



smt/i 
sin 6*' 
sm 62 



(A3) 



(A4) 




'")-!= cos ( 



, ^ smt/' ^ smw 

-2(ni + 1) cos — 2712 cos 6^2^ — 

sm di sm 62 



(A5) 



Looking at Fig. [19] it is clear that if cti < each excur- 
sion contains exactly one caustic, and the contribution /ij 
from the central part of the trajectory is exactly ni — ^2. 
If (Ti > 0, however, one of the excursions will be in the 
configuration schematized on the lower part of Fig. |19[ 
and will contain either two or zero caustics. This tran- 
sitional excursion can take place either in region 1 or in 
region 2. In the former case the number of caustics in 
the transitional excursion is zero if 61 < 7r/2 and two if 
9i > tt/2, and vice versa in the latter case. Since however 
caustics are counted with an opposite sign in regions 1 
and 2, this leaves the contribution unaffected. 

In the same way the number of caustics in the last 
excursion is one if <T2 > 0, and zero otherwise. Finally, a 
caustic can be found in the first excursion if Oi — 9' > tt. 
All this can be summarized by 

' m'J = e(a2) 

^} = m - 712 -t- sign(6li - 7r/2)e(CTi) (A6) 
^'j = e{9i -0' -tt) 
with Q the Heaviside step function. 

2. Transmitted trajectories 

This time, the final position is y" = R2 cos 9" — 
R2 cos ^2- Still working with y' fixed leads to the same ex- 
pressions as above for {dV /d9')y' and {dL" /d9')y: , while 
the variation of the final position can be shown to read 



-1 dy" 



89' 



sm ( 



(tan 61")- 



'(2771 + 1) cos 6*1 



y' 

sin ( 



sinf 



- (2712 + 1) cos 6*2 



cos ( 



smf 



(A7) 



leads to sn\9'd9' = sm9id9i, ai and a2 can be computed 



Defining the same indices as for the reflected trajecto- 
ries, the same expression is obtained for the number of 
caustics in the central part of the trajectory. Concern- 
ing the extremal portions, a caustic can be found in L' 
if 6'i — 6*' — TT > (as before) and one in L" this time if 
02 < 0. 
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Total Maslov index 



For semiclassical expressions such as Eq. (56), it is the 
sum /ij + lis which is relevant rather than fj,j alone. 
Although we will not provide a formal proof of this, it 
can be seen (and it is easily checked numerically) that 

f^j + f^s = - K''") + "1 - "2 + 1 • (A8) 
We make use of this equality in section |Vj 

Appendix B: Calculation of the prefactor J dy,S 

In this appendix, we compute the prefactor obtained 
in Eq. ( 56 ) once integral ( 43 ) has been evaluated in the 



stationary phase approximation. It reads 



12 



iy's) 



(Bl) 



Recall Sj is the action of the Green's function, Sn,a that 
of the mode, Jj the Green's function prefactor and y'^ 
the stationary phase point. Using basic properties of the 
action we have 

(B2) 

with x' = vf cos9n,a{ys) ^^nd .t" = vp cQsOn,a{y"), and 



Qyl2 



y" 



dy' 



(B3) 



y" 



Now, for a variation within the manifold on which the 
mode {n,a) is built, y" — y" {y' , p'y(y' , (n,a))), and thus 



P'y 



dp'y 



dy' 



dy' 



(B4) 



Using then that {dy" /dy')\ , 



{dp'yldy')\y. {dy"ldp'y)\ 



we obtain 



dy"dy' 



dy'' 



dy' 



(B5) 



which expresses the usual ratio between measures dy" 
and dy' on the corresponding Poincare sections. Making 



use of the identity W = {ni + l)Li + 712^2 + i?i(sin6'" 
sin 9'), it can be easily evaluated as 



= Ri smOn^aiy") 



dy" 


dy" 


de" 


de' 




dy' 
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„,a dy' 





COSen,a{y') \ ■ a ( 



and expression (|B1|) hence takes the final form 



,^ d^{Sj+S 

n.a J 
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(B6) 
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Appendix C: Main steps for the derivation of T„^a 

Here we describe how to deal with the case of trans- 
mitted trajectories (for which e(r") < 0) and obtain the 
transmission T„^a- We discuss only the few main steps 
of the calculation that differ from the derivation of -R„,q, 
presented in the main text. 

This time, the matrix structure of the integral that 
appears in Eq. ( 43 ) reads 




V-{r")V;^ir')a^ 



X cos 



(Cl) 



The Green's function projector V~ {r")Vj^ {r') selects 
those trajectories which connect the initial Poincare sec- 
tion in the positive (electron) eigenspace and the final 
Poincare section in the negative (hole) eigenspace. Ap- 
plying the stationary phase approximation to integral 
( 43 1 yields the same condition as the one obtained in the 



main text for the reflection, and the action of the Green's 
function reads Sj {y',) = hk^''' {x" - x') + (ni + 1 /2)S'^''" + 
{n-2. -I- 1/2) 5*2 v"S{y") -i^'SnAv's)- The Maslov index is 
computed as in the case of reflected trajectories (details 
can be found in appendix [A|. 



By putting together these elements, the integral (43) 
can be evaluated as 
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(5i"-°+S2)/2 



N2 



(C2) 



As expected (and contrary to the case of reflected trajec- 
tories), the incident mode once propagated to the point 
r" is no longer quantized. This, however, turns out to be 
irrelevant when it comes to computing the conductance 



making use of Eq. (42 ). Switching from position to angu- 



lar coordinates, the transmission probability of channel 
n polarized in valley a gives Eq. ( 62 ) . 
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